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Abstract This paper consists in two parts. First we set up a general scheme 
of local traps in an homogeneous deterministic quantum system. The current 
of particles caught by the trap is linked to the dynamical behaviour of the 
trap states. In this way, transport properties in an homogeneous system 
are related to spectral properties of a coherent dynamics. Next we apply 
the scheme to a system of Fermions in the one-particle approximation. We 
obtain in particular lower bounds for the dynamical entropy in terms of the 
current induced by the trap. 
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1 Introduction 

In this paper, we are interested in time scaling properties of propagation 
by coherent quantum dynamics. Non-trivial behaviour of a single-particle 
Hamiltonian can lead to a description of anomalous diffusion of electrons 
in solids [3]. This behaviour becomes apparent through the scaling of the 
spreading of one-particle wave functions with respect to time. We shall here 
adopt another approach: we introduce a localised trap in an infinite system 
and study the time behaviour of the current of particles falling in the trap. 
Applied to systems of Fermions in the non-interacting approximation, we 
obtain a lower bound on the dynamical entropy in terms of this current. 
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The trap states will be described by a collection of wave functions and we 
relate the current to the dynamics of these states. In particular, we show that 
an absolutely continuous spectrum produces a non-zero asymptotic current. 
A singular spectrum will lead to asymptotically vanishing currents possibly 
characterised by a dynamical exponent. 

A number of related topics and models have been considered in the liter- 
ature, mostly for the case of a continuous time evolution. The occurrence of 
singular continuous spectra as a source of anomalous diffusion is caused ei- 
ther by randomness in the Hamiltonian or by aperiodicity. We shall however 
not be concerned by producing such models but rather link scaling properties 
of dynamical entropy to assumed spectral properties of a discrete dynamics. 

Coherent transport in quantum systems is being studied by using reser- 
voirs as drivers. A number of delicate questions arise in this context with 
respect to the thermodynamic limit. Anomalous transport due to spatial 
randomness in the dynamics seems to occur [11]. 

Strongly chaotic classical or quantum dynamical systems generate en- 
tropy at a non-zero asymptotic rate: the dynamical entropy. In the classical 
case, the sum of the positive Lyapunov exponents is a bound for the entropy 
(Ruelle's inequality) and equality is reached for sufficiently smooth systems 
(Pesin's theorem). For quantum dynamical systems several entropies have 
been introduced such as the CNT construction based on a coupling with a 
classical system and the ALF construction that relies on POVM's (opera- 
tional partitions of unity). In order to obtain a non-zero entropy an abso- 
lutely continuous dynamical spectrum is needed, at least for Fermion systems 
in the one-particle approximation [10]. In open classical systems, the escape 
rate formalism links the escape exponent from an unstable repeller to dif- 
fusive transport. The escape rate is given by the missing exponents in the 
entropy for motion on the repeller [6]. 

There are however many mixing dynamical systems with less pronounced 
randomising properties which are not given in terms of exponents or rates. 
Such dynamics may lead to a sublinear scaling for the total dynamical en- 
tropy [4]. 

The aim of this work is to establish a lower bound for the entropy in 
terms of dynamical exponents of a localised trap in an infinite system both 
in the regular and the anomalous case. 

As a motivation we provide in Section 2 a few simple examples of the 
use of a trap in classical dynamics. Obviously there is a different physical 
mechanism at work with possibly similar macroscopic effects but our aim is 
to show that the time behaviour of the current at the edge of a localised trap 
encodes relevant information about the transport properties of the dynamics. 
Because of the locality of the trap, it is possible to deal immediately with an 
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infinite system and to avoid using a delicate large volume limit of boundary 
conditions. Section 3 deals with localised traps in unitary Hilbert space dy- 
namics and we study in Section 4 the dynamical entropy for non-interacting 
Fermions. We obtain in particular a lower bound for the entropy in terms of 
the current of particles falling in a trap. 



2 Traps in a classical context 

Consider first the simplest classical counterpart of a trap absorbing free elec- 
trons moving with velocities below that corresponding to a given Fermi level. 
This is a system of classical particles homogeneously distributed in space and 
with uniform velocity distribution below a maximal one. The state of such 
a system is a uniform mixture of spatially homogeneous states with fixed 
velocity. The number of particles per time unit caught in a localised trap is 
constant in time and essentially determined by the average cross-section of 
the trap. 

Next we consider the model of the diffusion equation in one dimension 
with a trap at the origin. We have to find the solution of the equation 

dn d 2 n 
dt dx 2 

for x > and t > with boundary condition n(0, t) — and initial condition 
n(x, 0) = 1. In this equation D is the diffusion constant and n(x, t) represents 
the particle density at time t and place x and the particle current is given 
by Fick's law 

d 

j(x,t) = -D—n(x,t). 
The solution of the diffusion equation reads 



>0 sin(kx) 
dk 



n(x,t) = ~ I - 
Therefore, the current at x = is 

^ Jo 



dke~ 



Dk 2 t 



-Dk'h 



Tit 



We recover hereby the usual exponent for diffusion. 

The third example is a simple random walk in one dimension with the site 
zero absorbing the walker. Let p(x, t) denote the probability that the walker 
reaches the origin for the first time at time t starting out at site x. We may 
assume that x G IN. The probabilities p are determined by the recursion 
relation 
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p(0, 0) = 1 and p(0, t) = for t > 



• p(x, t) — whenever x > t 

• p(x, = |(p(a: - 1, * - 1) + p(x + 1, t - 1)). 
It can be checked that the solution is given by 



p(x,t) = 




t — x even integer 
else. 



The current at time t is then given by 




In this formula, [x] denotes the largest integer less or equal to x. 

Finally, a random walk on Z 3 with a trapping set A leads to a total current 
Ja(£) = dNA(t)/dt where iVyi(t) is the total number of particles trapped by 
the set A up to time t. It turns out that 

J A (t) ~ C(A) + 2(27r)- 3/2 C(A) 2 t- 1/2 , 

where C(A) is the capacity of the set A [9]. Again, the behaviour of the 
current returns the relevant information on the transport properties of the 
system. 

3 Traps in Hilbert space 

We first consider an abstract model of a trap that absorbs particles. An 
explicit connection with a model of Fermions will be presented in Section 4. 
The basic ingredients of our description are an infinite dimensional Hilbert 
space $), a unitary operator U on Sj which specifies the evolution during a 
single time step, and a non-negative operator A less or equal than 1 which 
describes the effect of the trap. We shall assume that the trap is local in the 
sense that it is of finite rank d. Writing its eigenvalue decomposition 

d 

^ = I>^>(^-|, 0< Pj <l 

we can think of pj as the probability that a particle in the state \vpj) is 
absorbed by the trap when hitting the trap once. A value of pj close to 1 
means that the trap captures particles very efficiently in the state \ipj). 
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Starting out with a density < p < 1, the density after one time step 
and hitting once the trap becomes 

(1 - A)U pU*(l - A) = (TU)p(TU)* 

with T := 1 — A. Assuming a uniform initial distribution, i.e. p a scalar 
multiple of 1, the number of particles absorbed up to time t by the trap is 
proportional to 

N A (t) := Tr(j - . (1) 

The operator {TU) t {U*T) t is a finite rank perturbation of 1, therefore for- 
mula (1) makes sense. 

We shall always assume that eventually infinitely many particles are ab- 
sorbed by the trap. This will certainly be the case if the dynamics has 
reasonable randomising properties. We are in particular interested in the 
scaling behaviour of N A (t) with time, i.e., in the exponent 7 governing the 
asymptotics of N A 

N A (t) ~ t 1 " 7 , t large. 

The exponent 7 is non-negative as the growth of Na is at most linear in time 
and cannot exceed the value 1 by our assumption lim^ 00 A^^(t) = 00. We 
expect this scaling to be related to the spectral properties of U and to be 
independent of the size of the trap. Using the telescopic formula 

t-i 

1 - (TUY(U*TY = ^(Tf/) s (l - T 2 )(U*T) S , 

we rewrite for t > 1 N A (t) in terms of a current J A 

t 

N A (t) = J2Us) 

s=l 

with 

J A (t) := N A (t) - N A (t - 1) 

= Tr{ (1 - (TUY(U*TY) - (1 - (TUY^^TY' 1 ) } 

= Tr(Tf/)'- 1 (l-T 2 )(f/*T)*- 1 . (2) 

The expected scaling behaviour of J A is then 

Ja{^) ~ £~ 7 , t large. 
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The actual analysis will be performed for the simple case where A is a 
one-dimensional projector |y?)(y?| with tp a normalised vector in S) and we 
shall henceforth drop the subscripts of iV and J. In this case, N is fully 
determined by the probability measure 

IM(M) := \\E(d6M 2 (3) 

on the unit circle S l where E is the spectral measure of U 

U= f e- id E{d9). 
Js 1 

The expression (2) for the current now becomes 

J(t) = \\(l-P t - 1 )--.(1-P 1 )<p\\ 2 , t>l (4) 

with P t the orthogonal projector on U^ip. For consistency we must put 
J(l) = 1. As (1 — P t ) is a contraction, J is a monotonically decreasing 
function and therefore 

Joo := lim J(s) 

s^oo 

exists. The behaviour of the current provides information on the transport 
properties of the system. As we don't have in our general description a notion 
of position operator, allowing a definition of ballistic or diffusive motion in 
terms of spatial spreadings of wave functions, we shall rather concentrate 
on the relation between the current and the randomising properties of the 
dynamics which are quantified by the dynamical entropy. This is, within 
the context of non-interacting Fermion systems, the subject of Section 4. In 
this section we shall study the asymptotic current in terms of the dynamical 
properties of the trap states for the simple one-state trap. The analysis could 
be extended to more involved traps that have a spatial structure. The more 
complicated behaviour of the subsequent current could then be studied as 
in the case of the 3D classical random walk, leading possibly to quantum 
capacities. 

In order to relate the measure \i with J^, we introduce for \z\ < 1 the 
function 

oo 

(?(*):= 5>V(«), (5) 

s=l 

where /i A is the Fourier transform of \i 

^(t) : = I f x(d9)e- lW , t E Z. (6) 
Js 1 



6 



Expressing G in terms of /x, we find 



Obviously, G is analytic in the open unit disc and we shall be concerned with 
its value on the boundary of the disc. Writing z = re tv with < r < 1, a 
direct computation shows that 

1 + 2BeG(re i ") = / /x(d0) — 2 1 ~ - ^. (7) 

y^i 1 + — 2r cos(?] — 0) 

The function 

1 -r 2 

* r (0) : 



1 + r 2 — 2r cos6 7 

is the Poisson kernel and the S r are a 5-convergent sequence of smooth, 
positive, normalised functions 



J /-27T 



2tt 

and 



d#<5 r (0) 



/fo) =limi / Mf(e)6 r (6- V ) a.e. (8) 

for any integrable function / on the unit circle. By a.e. we shall always mean 
almost everywhere with respect to the Lebesgue measure. 
The imaginary part of G is given by 

SmG(re-) = I /x(d0) ^ " &) 



l s i 1 + r 2 — 2r cos(r/ — 9) ' 

When r tends to 1, we obtain the Hilbert transform of /x [5, 7] 

sm(i] — 9) 
1 + r 2 — 2r cos(i] — 9) 



(W/i)fa)=liin/ /x(^) 



= limi/ /x(d0) cot( V^)- (9) 
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The limits in (9) exist almost everywhere and for each e > the set on which 
|7Y/x| is larger than 1/e has Lebesgue measure less or equal to e. We shall 
now express the current and asymptotic current in terms of G and thus in 
terms of the spectral properties of the trap, i.e. of the measure /x. 
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Lemma 3.1 With the notation of above 



t-i 



J(t) = l-J2\K(s)\ 2 , 



s=l 

where the function K is determined by the relation 

v ' s=l 

The asymptotic current is given by 



'oo 

rT 



Ti 2vr Jgi |1 + G(re 



i8\\2 ' 



Proof: 

We introduce 



K{t) :=((U*) t y>,(l-P t _ 1 )---(l-P 1 )<p) 
and denote by F the Z-transform of K 

oo 

F(z) := £Vtf(s), \ z \< 1 - 

8 = 1 

The relation (10) follows from a straightforward computation 

j(t) = ||(i -p t _o ...(i- po ^|| 2 

= iki - p t _ 2 ) . . . (i - po v^r - K(c/*rv, (i - Pt-2) • • • (i - Pi 
= j(t-i)-i^-i)i 2 . 

Next, we write 

(I _ Pt ) . . . (J _ Pl ) y = (J _ p t-1 ) ...(!_ Pl ) ^ _ # (j) ([/*)V. 
Therefore 

t 

(1-P t )...(l-P 1 )<p = <p-Y I K(8) (tT)V- 

Taking the scalar product of (£/*)' +1 y? with (14), we obtain 

t 

Kit + 1) = ((frrv v> - £<(tn^ + v, *>> *w 

s=l 

8 



Equation (15) has the structure of a one-sided convolution equation. Using 
the Z-transform it becomes an algebraic equation. More precisely, multiply- 
ing (15) with z t+1 and summing from t = to oo, we obtain 

F(z) = G(z)-F(z)G(z). 

As by (7) 1 + G never vanishes inside the unit disk 

J^L. F(z) = ±, K(s) , 

v ' s=l 

proving (11). 

The basic relation between the asymptotic current Jqo and the measure \i is 
obtained by applying Parseval's formula to (13). Putting z = re ir > 



1 

2^ 



S 1 



Fire 



J2r 2s \K(s)f. (16) 



s=l 



We compute now the asymptotic current on the basis of (10) 



t-i 



Joo = lim J(t) = 1 - lim V \K(s) 

t^oo t^oo ' * 

s=l 



oo „ 

= l-limVr 2s |K(s)| 2 = l-lim— / d6\F(re t6 )\ 2 (17) 

r|l ^— ^ r|l 27T J51 



1 f l+2KeG(re ie ) 
= lim — / dB 



rTi 2tt J s i |1 + (^(re* 1 
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Our first result deals with the asymptotic current for trap states with 
absolutely continuous dynamical spectrum. 

Theorem 3.1 Suppose that /i is absolutely continuous w.r.t. the Lebesgue 
measure, then > 0. 

Proof: 

Let jj, be absolutely continuous with respect to the Lebesgue measure with 
density p. We have for r < 1 

(i + 2 S ec) (re »o = jl j Mp( e + ,) 1 _ 2 ^ ro ;; +r2 . 
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Because p is integrable and because of the properties of the Poisson kernel 

lim(l + 2M)(re") = p(rj) a.e. 

rfl 

Also the imaginary part of 

z^ — [ d9p(9)-^— 
converges almost everywhere to the Hilbert transform Hp of p(6)d9: 
limSmG^re'") = (Hp)(r}) := lim / d9 p{9 + rj) cot(0/2). 

rTl SIO 47T J\ d \>g 

In the expression for the asymptotic current 
Joo^lim^- / d9(l - \F(re i0 ^ 



s 1 



1 f ,„ l + 2&G(re f 

hm — / d0 ■ 



rti 2tt y 5 i 1 + 2KeG(re ie ) + (3?cG) 2 (re lf? ) + (9fmG) 2 (re <e ) 
the integrand is bounded by 1 and tends almost everywhere to 

(l + p(0))2 + 4(ftp) 2 (0) 

as r grows to 1. We can therefore apply the dominated convergence theorem 
to obtain 

t 1 f AO , n 

' ~ 2n] s ™ {i + p{9)y + A{H P y{9) >u - 



In (17), we have replaced a limit £ — > oo by a limit rfl. In fact more 
information can be gained in doing so. Indeed, in the case Joo = 0, the 
behaviour of J(£) for large £ is that of 1 — X)!=i l-^( s )| 2 - The function 

oo 

J{r):=l-Y j r 2s \K{s)\\ 0<r<l (18) 

is expressed in terms of the discrete Laplace transform of s \- > |i^(s)| 2 and 
Tauberian theorems relate the large time behaviour of £ i— > J(£) = 1 — 
Xls=i l-^( s )| 2 w ^h that of J(r) when rfl. We shall first use this idea to 
show that Jqo vanishes when p is singular. Next, we shall in a few examples 
consider convergence exponents. 
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Theorem 3.2 If ji is singular w.r.t. the Lebesgue measure, then Joo = 0. 



Let /i be singular and hence concentrated on a measurable set with zero 
Lebesgue measure. As the Lebesgue measure of this set equals the infimum 
of the Lebesgue measures of open subsets containing the set, we can given any 
positive e find an open subset A of S 1 such that the Lebesgue measure of A 
is not larger than e and fi(A) = 1. The set A is a countable union of disjoint 
open intervals ]aij, f3j[. We dress each of the ]aij, f3j[ with open strips of width 
5j which shall be determined later on. By choosing the 5j sufficiently small 
we may still ensure that the Lebesgue measure of A + := [_}j]oij — Sj,f3j + Sj[ 
is small. 

We shall also need 



Proof: 




for 5 < 1 — r 



1 - 
< 

- 5 



for S > 1 — r. 



(19) 



The estimate for the case 5 > 1 



r is obtained as follows: 




\e\>s 



I _|_ r 2 _ 2r cos 




I _|_ r 2 _ 2r cos 




2(1 -r) 




1 



7T 



(l-r) 2 + 2r(l-cos#) 





for S sufficiently small. 
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We now estimate the current 



A r ) = ^ J si d V (l ~ \F(re">)\ 2 ) 

< e + 2^ 5 , + ^/^ + d,(l-|^)| 2 ) 

- € + 2 ^ + 2n J S1 \ A+ dV 1 + 2KeG + (KeG) 2 + (3mG) 2 (re } 



l + 2KeG , in , 
' re" 



" 6 + 2 ^ ^ + 27 1 w ^ T + 2KeG + (^eG) 2 
< e + 2V 6j + - [ d V [ n{69) - ' 

^ 3 7t j sl , A+ y 51 PV i + r 



/si\a+ Js 1 ' 1 + r 2 - 2r cos(r? - 9) 

< 6 + 2Y 5j + -[ pi(d9) [ dr}S r (rj-9). (20) 

In the last but one inequality we have used that the integrand is bounded 
from above by 4(1 + 29?eG). In the last inequality we use (19) to get 

J(r)<e + 2 J £5j + l( E MW^D) 

j k 

S k <l-r 

+ 2 (1 _ r) / y KHM\ (21) 

7T V ^ 6 k J 

& k >l-r 

For a given e and a given large N, we may take 5k = N/iQotk, — ?")■ 

When r is sufficiently close to 1, the upper bound for J becomes 

J(r)<e + 2iV(l-r) + ^( £ MH,A[)) 

A r M(]"fe fe &[)<l 

+ ^#({* WK,ftD>i})- (22) 

First we fix an arbitrary small e and the corresponding sets ]ccfc,/3fc[. The 
last two terms in (22) can be made small by choosing iV sufficiently large. 
Finally the second term in (22) becomes small when we let r f 1. 
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We conclude this section with some examples, remarks and partial results 
about exponents. Let us assume that we are in the situation = and 
that we can assign an exponent to J, i.e. 

7 := , im _l!*£W 

t^oo log t 

exists. We shall, moreover, assume that infinitely many particles eventually 
are absorbed by the trap and actually strengthen this condition to 7 < 1. 
We are interested in relating the behaviour of t 1 — > J(t) as t — > 00 with that 
of r 1— > J(r) as r | 1. We therefore introduce upper and lower exponents a 
and a for J 

_ log J(r) . log J(r) 

a := limsup - — and a : = limmt 



rT i log(l-r) _ ' r|i log(l-r)' 

Lemma 3.2 [/sing £/ie notations and assumptions of above, a = a — 7. 
Proof: 

We have to consider 

t— 1 00 
J(t) = l-^|K(s)| 2 = ^|^(s)| 2 and 

s=l s=i 

00 00 

J(r) = 1 - \K(s)\ 2 = " r 2 W*)| 2 - 



s=l s=l 



For notational convenience, we put c t : = \K(t)\ 2 and A := — 21ogr. 
Let a > 7, then 



00 1 00 

s=l ' ' s=[A-!] 

00 

= lim (1 - -) N a V c s = 



OO. 

N— >oo" e' z — ' 

s=JV 

Hence, a > 7. 

Conversely, let ct < 7 and introduce the short notation f(N) := J2^L N c s . 
Fix an arbitrary e > and a 70 such that a < 70 < 7. For A > and 
k = 2, 3, . . ., let N n be determined by 

i \ i\~t \ k — n 
exp(— \N n ) = , n — 1, 2, . . . , k — 1 
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and put N n := [N n ] + where [x} + is the smallest integer larger or equal to 
x. We shall pick k later on in such a way that k <C A -1 . This implies that 
the Nj are far apart, in particular that Ni >> 1. We now have obtained a 
partition 

1 =: iVo < ^ < • • • < iV K „! < JV K := +00 

of 1N such that 

n 4- 1 

1 _ e ~ A * < — 1_ for iV n < t < N n+i . 

K 

OO K — 1 iV n + l — 1 

D 1 - e ">- = E E (!-e- As ) Cs 

s=l n=0 s=AT n 

, 1 N n+1 -1 

E - 

n=0 s=N„ 

= E— (/W-Z^n+l)) 
z — ' ft 

n=0 

<-(/(iv ) + /(iv 1 ) + ... + /(iv K _ 1 ; 

< - f 1 + eTVr 70 + • • • + eAT 7 ? } . (23) 

The inequality in (23) holds for N± large enough, i.e. for k sufficiently large. 
Because N n = log(^^-)] + we have 

JV^ < - log ' 



A V 1 - 2 



and thus 



E(1 _ e -,. )c ,<I + I £E (l log (^)) 

s=l n=l v v « / / 

n=l K 

< - + A 7 °e [ dy flog(- 



-70 



= - + A 7o e / di flog-) 

K J \ tJ 

= - + A^er(l- 7o ) 
< 5A 70 



1 \ _ T0 
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with 5 arbitrarily small. This last inequality is obtained by choosing A~ 70 <C 
k <C A" 1 which is possible because < 70 < 1. Hence a < 7 and the lemma 
is proven. 

■ 

When estimating the current in the proof of Theorem 3.2 we dropped 
the contribution of SmG. Generally, this may lead to underestimate the 
exponent. A simple example is provided by a measure \x that is concentrated 
on a finite set such as p(d&) = 5(6)d6 (the general case being quite similar). 
A simple calculation shows that F(z) = z. Therefore 

— / dr] (l - \F(re tri )\ 2 ) = 1 - r 2 ~ 2(1 - r) 
2n J s i \ J 

and the true exponent is 1. Dropping the imaginary part of G in the integral, 
we get an exponent 1/2: 



2tt J s i 1 1 + 2KeG + (KeG) 2V ' 2ix J s i 1 (l-rcos^) 2 

[l-r)x 2 



I 



dx 



(2(l-r)+x 2 ) 2 
~ y/l — r. 

Let \i be a discrete measure, possibly concentrated on a dense subset of 

fj,(de) = J2pj 6 ( 6 - e M e 

3 

with 

Pj > 0, pj = 1 and i ^ j ^ di ^ 6j. 

j 

Applying the estimates in the proof of Theorem 3.2 we obtain 

J»<25> + 4 ^ ^=^ + 4 £ (24) 

<5j>l-r 5j<l-r 

Suppose that the pj tend sufficiently rapidly to zero in order that also 
J2i y/Pj < 00. Choosing 5j = a/(1 — r)pj, we obtain J(r) < \fPlj) 



j y rj - ~~- ~ — — o "j v v~ /fJ> "~ " v / — V rjj v - r 

and therefore an exponent 1/2. In view of the first example, this is the best 
exponent we can hope to obtain neglecting the contribution of the imaginary 
part of G. 
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Suppose that pj ~ j a with a > 1. Choosing in (24) 5j = Npj(l — r) 
with N large, we obtain 



The optimal choice for iV is N ~ (1 — r)~ a ^ 2a ~ 1 ^ and this yields an exponent 
(a — I) /(2a — 1) < 1/2. The exponent 1/2 is reached for all pj that tend 
faster to zero that any inverse power of j. 

In our last example the measure p will be singular continuous [8]. Given 
a number x G [0, 1], we write its binary expansion as 



where a m (x) G {0,1} for m > 1. This defines a map F : {0,1}™ -> [0,1] 
which can be used to transport a measure A on {0, 1}^ to a measure p on 
[0, 1] by the relation p(A) = X(F~ 1 [A]). Take for A the infinite product of 
the measure (1 — p,p), p G [0, 1]. The corresponding measure p on [0, 1] is 
called the Bernoulli measure p p . Except for p = 0, p = 1 (Dirac measures) 
and p — 1/2 (Lebesgue measure), p p is singular continuous. 

We want to estimate the exponent a of the current J(r) ~ (1 — r) a . 
Assume < p < 1/2 and let q be a number between p and 1/2. Define the 
sets A(n) for n > 1 as 

A(n) = {x G [0, 1] | #{m < n \ a m (x) = 1} < gn}. 

The asymptotic behaviour of the Bernoulli measure is 

1 - p p (A(n)) ~ exp(-nS(q\p)), 

whereas for the Lebesgue measure 

|A(n)|~exp(-nS((z|l/2)), 

with S(pi\p2) the relative entropy of the probability measures (pi, 1 — pi) 
with respect to (p2, 1 — P2), i-e., 

S(pi\P2) = Pi log(pi) + (1 - pi) log(l - pi) - pi log(p 2 ) - (1 - pi) log(l -p 2 ). 

The sets A(n) are finite unions of K(n) intervals A k (n) for which the first n 
digits of the binary expansion are given. We have the asymptotics 



J(r) < 2N(l - r) + ^{{j I Pj > ^}) + \ E ^ 




00 



K(n) ~ exp(nS'(g)) 
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where S(q) is the (Shannon) entropy of the probability measure (q, 1 — q), 
i.e., 

S(q) ■= -glog(g) - (1 - q) log(l - q). 

As in the proof of Theorem 3.2 we dress the intervals A k {n) by small 
strips of length 5k{n). This results in the set A + (n). We can then use the 
estimate (20) for the current J(r) 



dry I [ LL(d9)5 r (ri -9)+ [ v(d8)5 r (r} - 6) 

{jA(n) JS\A(n) 



S\A+(n) KJA(n) J S\A(n) 

<- M^(n)) + (l-r) + 

5 k (n)<l-r 5 k (n)>l-r 

Again taking Sk(n) = Nfj,(Ak(n))(l —r), we obtain 

J(r)<\A(n)\+2N(l-r)+ ^ fi(A k (n)) 

k 

N^(A k (n))<l 

+ 1#{A; | N^(A k (n)) > 1} + (1 - pi p (A(n))). 

Imposing now the scaling behaviour N ~ exp(nz/) and 1 — r ~ exp(— n/3), 
we find 

- mm{S(g\l/2),P-u,u-S(q),S(q\p)} \ 

for any q, v and (5. The optimal values are 

log 2(1 -p) 



a > 



Q = 



° l— p 



(for which S(q\l/2) = S(q\p)), v = log 2 and (3 = 2 log 2 - S(q). Finally, the 
lower bound for the exponent a is 

log 2 -5(g) 

" £ 21o g 2-S( 9 )' (25) 

We compare this bound with some numerical computations. For a few 
pairs (p, r) the integrals for G(z) and J(r) were evaluated on an equidistant 
mesh of 2 13 points. To illustrate the approximation made in (20) by dropping 
$5mG, we performed the computation both with and without this imaginary 
part. The relative precision for the current J(r) was checked to be better 
than 0.01, while for the exponent a it is 0.1. Figure 1 shows the existence of 
the exponents and the error introduced by neglecting the imaginary part of 
G(z). The latter is quantified in Table 1. For this example the lower bound 
(25) for a is rather sharp. 
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Figure 1: The current J(r), with (crosses) and without (circles) the contri- 
bution SraG. Left, p = 1/3 and right, p = 0.95. 





analytical 


numerical 


numerical 






without ^smG 


with SmG 


p=l/3 


2.05 10~ 2 


3.7 10- 2 


5.6 10- 2 


p = 0.95 


1.96 10- 1 


3.2 10- 1 


4.2 10" 1 



Table 1: Estimates for the exponent a. 
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4 Dynamical entropy of a Fermion dynamics 



We shall in this section apply our results in the setting of a free Fermionic 
gas. As this is a system of non-interacting particles, it is completely de- 
scribed in terms of single-particle quantities. Second quantisation allows to 
lift one-particle objects to the many particles, taking into account the Fermi 
statistics. We remind here briefly the mathematical setup. 

We shall denote the single-particle Hilbert space by Sj. The observables 
of the Fermion algebra, also called CAR for canonical anticommutation re- 
lations, is the C*-algebra A($)) determined through the relations 

a(f + ag) = a(f) + aa(g), {a(f),a(g)} = and {a(f), a*(g)} = (/, g). 

Sometimes we shall deal with a one-particle space f) of finite dimension d. 
In this case A(Sj) is easily seen to be isomorphic to the algebra of matrices 
of dimension 2 d . An explicit construction is given in terms of linear trans- 
formations of the antisymmetric Fock space T(Sj) which is spanned by the 
n-particle vectors 

a*(h)a*(f 2 )---a*(f n )Q 

for < n < d. The normalised vector Q is called the vacuum and it is 
annihilated by any operator a(f). 

The construction of dynamical entropy presented in [2, 1] is based on the 
following idea. Given a unital C*-algebra 21 and a reference state u, one 
considers an operational partition, i.e. a finite collection X = {xi, x 2 , ■ ■ ■ , x n } 
of elements of 21 satisfying Ylj x j x j = !■ This yields a correlation matrix 

p[X] := [u(x*Xi)\ 

with corresponding von Neumann entropy 

H[u,X] :=Tr^(p[X]), 

where r\ is the usual entropy function 

i](0) := and r](x) := — xlogx, < x < 1. 

The average entropy of an operational partition X = (x±, x 2 , ■ ■ ■ , x n ) arises by 
computing the average entropy of the correlation matrix % t corresponding to 
the refinement of the partition 3C at discrete times up to t — 1. More precisely 

X t := e W (I)o...o9(I)ol 
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In this expression, 



S (X) := (e s ( Xl ), s (x 2 ), . . . , e s (x n )) and 

X o 2) := (xij/i, x 2 yi, . . • , x„y m ) for 2) = (y u y 2 , . . . , y m ). 



It can happen that the growth of t i— > H[c<j,£i] is sublinear if the dynamics 
is not sufficiently randomising. In such a case one can look for a growth 
exponent. 

We shall in the following pages obtain a lower bound for H[c<j, £ t ] in terms 
of particle numbers absorbed by a trap for the case of a weakly interacting 
Fermion system, meaning that we may use an effective one-particle dynamics 
6(a(/)) := a(Uf) for the evolution. U is a unitary operator on the one- 
particle space Sj. 

The reference state will be chosen accordingly as a gauge-invariant quasi- 
free state. Such a state ojq is uniquely determined by its symbol Q which is 
a linear operator on satisfying < Q < I. The only monomials in the cre- 
ation and annihilation operators a* and a which have non-zero expectations 
contain a same number of each and 



In particular, we may choose Q = k1 for < k < 1. Such states are 
homogeneous and in e.g. the case of Fermions on a lattice they describe 
independent Fermions occupying each site of the lattice with probability 
k. The Fock vacuum, i.e. the vector state determined by VL in the Fock 
representation of above, corresponds to the choice k — 0. The choice k—1/2 
corresponds to the unique tracial state on A(Sj). A quasi- free state ojq is 
known to be pure if and only if Q is an orthogonal projector. Moreover, 
any uoq can be obtained as the restriction of a pure quasi-free state on a 
larger CAR algebra by using the purification construction. One introduces 
the auxiliary space R :— Q(l — Q)Sj and the projection operator 



on S) © R. For the homogeneous states Q = k1 with < k < 1, & = $) and 
the projector becomes 



^<j(a*(/i) ' ■■ a *Un)a(g n ) ■ ••a{g 1 )) = det((&, Q fj)). 




(26) 




) 
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For a symbol Q of finite rank, the entropy of ujq is given by 

S(u Q ) = Tr(ri(Q)+ri(l-Qj). (27) 

The formula can obviously be extended to compact Q with eigenvalues con- 
verging sufficiently fast to 0. 

In order to compute the dynamical entropy for quasi-free evolutions with a 
quasi-free reference state, it suffices to consider a restricted class of partitions 
X characterised by the property that 

V i-> A(y) :=^2x*jyxj 
j 

transforms the gauge-invariant quasi-free states into themselves. Such maps 
A are called gauge-invariant quasi-free completely positive maps and are de- 
termined by two linear operators V and W on obeying the restrictions 

< W < 1 - V*V. 

On a monomial A acts as 

A(a*(f 1 )---a*(f n ))= e(S)(j[a*(Vf j ))u; w ( K l[a*(f k )). 

Sc{l,...,n} jeS k£S 

In this formula, a* denotes either a or a* and e(S) equals ±1 according to the 
parity of the permutation defined by S. The quasi-free state Uq transforms 
under A into the quasi-free state with symbol 

V*Q V + W. (28) 

Even if different partitions may yield the same map A, H[u;q,£] will only 
depend on Q and A and we shall derive in the next proposition its expression 
directly in terms of Q and A. 

Proposition 4.1 Let dim(^) < oo and let X = (xi,X2, ■ ■ ■ ,x n ) be an oper- 
ational partition in A(Sj) such that y i— > Ylj x *jy x j ^ s gauge-invariant quasi- 
free determined by (V,W). Let the symbol Q determine the gauge-invariant 
quasi-free state ujq, then 

H[u Q ,X] = S(co R ) (29) 
where R is the symbol on S) © (Q(l — Q)Sj) =: $) © A given by 



R 



V*QV + W V*y/Q(l-Q) 
VQ(1-Q)V (1-Q) 



si 



(30) 
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Proof: 



Let us denote by (S)q, tiq, Qq) the GNS triple of ujq and by (ei, e 2 , . . . , e n ) 
the canonical orthonormal basis of C n . The pure state on 23 (fig) <E> -M n 
induced by the vector Ylj^Q^^^Q ® e j £ £)q ® C n restricts to generally 
mixed states on M. n and *8(.?)q) that have, up to multiplicities of 0, the same 
spectrum. A straightforward computation shows that this restriction to M. n 
is the correlation matrix p[X] and that to *8(.?)q) the density matrix 



As A(Sj) is isomorphic to the algebra A / i 2 dim (-s) we can write that Sjq = C 2 "®£ 
with ttq(x) = x <g> l£. Therefore, there exists a unique unity preserving 
completely positive map T on *B(5]q) determined by the requirement 



Obviously H[co»q,£] = S(\Qq)(Qq\ o T). It now remains to compute this 
quantity. 

Using the purification (26) we see that the GNS representation space of ujq 
is the Fock space built on S) © (Q( J — Q)Sj) and that T is determined by the 
operators (V © 1, W © 0). It suffices now to use formulas (28) and (27) to 
finish the proof. 

■ 

We shall now obtain a lower bound on the dynamical entropy in terms 
of currents of particles falling into a trap. This will generally not provide 
the optimal lower bound but we expect it to provide the correct growth 
exponent, which it certainly does in the case of linear growth. The second 
quantised version of a localised trap is provided by a quasi-free completely 
positive map A with operators (V,0). In order to compute the number A 
of particles that disappear from an homogeneous state u K in the trap, we 
consider the particle number operator N in A(S)). N := Xlj a *( e i)°( e i) 
where {e 1 , e 2 , . . .} is an orthonormal basis of ft. We assume for the moment 
that $) is finite dimensional but the general case can be obtained by a suitable 
limiting procedure. Then 



M x j)^Q)M x j)^Q\- 



j 




A = cu K (N-A(N)) 



k Tr( J — V*V). 



The locality of the trap is expressed by the condition 



Rank(l - V*V) < oo, 



V*V < I. 
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In our case, all refined and evolved partitions remain quasi-free with 
strictly local action. An explicit computation shows that X t is determined 
by (V t , 0) with 

V t = [V [/]*[/-*. 

Using the explicit expression of the entropy of correlation matrix (29) in 
terms of its symbol (30), we have 

H[u K ,X t ]=T±(r j (R t )+r j (l-R t j) (31) 

with 

/ Kv;v t v^T^Jy; \ 

Rt -{^0^)v t 1-k )■ (32) 

In order to avoid a trivial situation we assume that < k < 1, in which 
case the trace in (31) is taken over a space of dimension twice the rank of 
1 - V t *V t . We write 

KV*V t ^k(\ - k) V t * 
a/ k(1 — k) V t 1 — K 



v / r^ o ) 1 o o 



A* A. 



But then 



A A* 



l-K + KV t V* 




Using that, up to multiplicities of 0, A* A and A A* have the same spectrum 
and that 77(0) = r](l) = 0, the entropy becomes 

Hk, **] = ttt7(i - «(i - u;u t )) + 77(«(i - u;\4)). 

Finally, as 77 is concave we obtain the lower bound 
Proposition 4.2 

Hk.lt] > MK)+ri(l-K)}Tr(l-V t *V t ). 

The usual computation of dynamical entropy involves two more steps. 
First the computation of the asymptotic rate of entropy production h[u;,X] 
which consists in taking the limit for t — > 00 of Hk,Xj/t and next taking 
the supremum over a suitable class of operational partitions. Proposition 4.2 
is more general in the sense that it provides a lower bound for H[o;,Xt] even 
when this quantity scales in a sublinear way in t. We can also use the 
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lower bound of the proposition to show that the dynamical entropy is strictly 
positive whenever there is a non-zero asymptotic current for a trap belonging 
to the class of allowed partitions. We have seen that k Tr(I— V t *V t ) represents 
the total amount of particles that disappeared from the homogeneous state 
uj k into the trap up to time t. The corresponding current at time t is then 

K-nwvt-vuvt-il 

which is precisely the quantity considered in Section 3. In particular, the 
entropy grows linearly in time if the absolutely continuous spectral subspace 
of the single-step unitary U is non-trivial. But even if U has no absolutely 
continuous spectral component, an estimate of the growth exponent of the 
entropy may be obtained. 
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